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1 Introduction 

The three-gluon coupling is perhaps the most obvious manifestation of the non-Abehan 
aspect of Quantum Chromodynamics (see also the reviews @, 0). Implicitly, it has 
been studied experimentally through the observed running of the coupling constant [Q. 
The associated Casimir invariant has even been measured directly in studies of four-jet 
events at LEP the SU(3) group being consistent with the data. 

Apart from being a standard object of consideration in textbooks on Quantum Field 
Theory and QCD (see e.g. P, ^ §), the perturbative corrections to gluonic vertices 
are also very important in real physical calculations, such as multijet production at the 
hadron colliders (see e.g. |P, |^ and references therein). At the present level of accuracy, 
one needs to perform not only calculations with on-shell external particles, there are also 
contributions where general off-shell results are needed. 

One of the original reasons the three-gluon vertex was studied was the belief that its 
infrared properties might shed light on the mechanism of confinement. In these stud- 
ies, different approaches were used, some of which are discussed in the review flOl (and 
references therein). 

For special cases, the one-loop results for the three-gluon coupling have been known 
for many years. Celmaster and Gonsalves (CG) presented in 1979 the one-loop result 
for the vertex, for off-shell gluons, restricted to the symmetric case, pi = p\ = p\i in 
an arbitrary covariant gaugeQ. Ball and Chiu (BC) then in 1980 considered the general 
off-shell case, but restricted to the Feynman gauge |jl3|. Later, various on-shell results 
have also been given, by Brandt and Frenkel (BF) |]T^ , restricted to the infrared- singular 
parts only (in an arbitrary covariant gauge), and by Nowak, Praszalowicz and Slomihski 
(NPS) [jl5l, who also gave the finite parts for the case of two gluons being on-shell (in 



Feynman gauge). An overview of these results is given in the table below. 





all momenta off-shell 


some momenta on-shell 


general case 


2 2 2 

Pi=P2= Pi 


pi = 


pI=pI = {} 


Feynman 
gauge 


BC0, 
eq. (3.3), 
no quarks 


special case 
of CG [0] 


special case 
of BF |14| 


NPS p, 
App. B 


Arbitrary 

covariant 

gauge 




CG ii, 

eq. (14) 


BF 0, 
eq. (25), 
no quarks, 
no finite parts 


BF in, 
eq. (30), 
no quarks, 
no finite parts 



Table 1: Kinematics and gauges considered in other studies 



From this table one can see that, even if we consider the results in (or around) four 
dimensions, there are still several "white spots". They correspond not only to the most 

"'The result of was also confirmed by Pascual and Tarrach |l2j . 
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general case (the lower left corner), but also to some other cases when the results are miss- 
ing, either for quark loop contributions or for the finite parts. The aim of the present paper 
is to cover all such remaining spots (for the case when massless quarks are considered). 
Moreover, we present results which are valid for an arbitrary value of the space-time 
dimension. Apart from the three-gluon vertex itself, we also consider the ghost-gluon 
vertex and two-point functions, to be able to check that all these expressions obey the 
Ward-Slavnov-Taylor identity for the three-gluon vertex. 

At the one-loop level, the simple and well-known Lorentz structure of the lowest-order 
coupling gets modified. In the general case, six tensor structures (and their permutations) 



are needed to decompose the three-gluon vertex ||T3[. Thus, six scalar functions multi- 
plying these tensor structures are to be calculated. These scalar functions depend on the 
gauge parameter, the space-time dimension, and the kinematical invariants {pf, p\, p^. 

There are several reasons why the one-loop results calculated in arbitrary gauge and 
dimension n are of special interest: 

(i) knowing the results in arbitrary gauge, one can explicitly keep track of gauge invariance 
for physical quantities; 

(ii) if one is interested in the two-loop calculation of the three-gluon coupling, one should 
know one-loop contributions in more detail; 

(iii) results in arbitrary dimension make it possible to consider all on-shell limits (when 
some pj = 0) directly from these expressions (see section 4), this is impossible if one only 
has the results valid around four dimensions; 



(iv) QCD is also a theory of interest in three and two dimensions (see e.g. []T6[ and the 
review [p!7[| ); 

(v) as we shall see, the results for arbitrary dimension are not much more cumbersome 
than those considered around four dimensions (in some respects, they are even more 
transparent and instructive). 

We note that in several papers the one-loop three-gluon vertex in axial-type gauges 
(including the light-cone gauge) was considered (mainly divergent parts and special 



limits have been studied). The three-gluon vertex in the background field formalism was 



considered in ref. W^, while the gauge-invariant vertex was studied in ref. |20|. Moreover 



there were some lattice calculations of the three-gluon vertex, see e.g. ref. |PT[. We shall not 



address these issues here, but concentrate instead on the standard vertex in an arbitrary 
covariant gauge. 

The paper is organized as follows. In section 2, we introduce the notation for the two- 
and three-point functions to be considered, and discuss their decomposition in terms of 
scalar functions as well as the corresponding Ward-Slavnov-Taylor identity. In section 3, 
we present the most general off-shell results for the three-gluon vertex. Section 4 con- 
tains the corresponding expressions for all on-shell limits of interest. In section 5, we 
conclude with a summary and a discussion of the results. Then, we have several appen- 
dices where some further results and technical details are presented, such as the formulae 
used to decompose the three-gluon vertex (Appendix A), relevant results for the scalar 
integrals involved (Appendix B), results for the self energies (Appendix C) and the ghost- 
gluon vertex (Appendix D), expressions for the on-shell limit p| = in arbitrary gauge 
(Appendix E), and also some results for Pi= P2 = ^ (Appendix F). 
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2 Preliminaries 

The Yang-Mills term of the QCD Lagrangian yields the following well-known expression 
for the lowest-order three-gluon vertex: 

- i 9 P'"'"'^ [^MiM2 (Pi - P2)m3 + ^M2w(P2 - Pa)/.! + ^A.3Mi (P3 - Pi)^2] . (2-1) 

where pi, P2 and are the momenta of the gluons, all of which are ingoing, P1+P2+P3 = 0. 
In ( p.l|) , the /"i"2a3 Qj-Q totally antisymmetric colour structures corresponding to the 
adjoint representation of the gauge group^. They can be extracted from the general 
three-gluon vertex by defining 

r™3(pi'P2,P3) = -i g r^'^' r,,,,,,{p,,p,,p,). (2.2) 

Since the gluons are bosons, and since the colour structures /'^i"2a3 qj-q antisymmetric, 
^fiiij.2fj.3iPijP2,P3) must also be antisymmetric under any interchange of a pair of gluon 
momenta and the corresponding Lorentz indices. 
The lowest-order gluon propagator is 

^"""^(w^-e^). (2.3) 

where ^ is the gauge parameter corresponding to a general covariant gauge, defined such 
that ^ = is the Feynman gauge. Here and henceforth, a causal prescription is under- 
stood, + iO). 

When one calculates radiative corrections to the three-gluon vertex (the corresponding 
one-loop diagrams are presented in fig. 1), other tensor structures arise, in addition to the 
lowest-order expression ( |2.1|) , and the general tensor decomposition should be considered. 
If we take into account momentum conservation (only two of the external momenta are 
independent), fourteen independent tensor structures carrying three Lorentz indices exist, 
and in general r^^^2At3 "^^^ t)e written as a sum of these tensors multiplied by scalar 
functions (see Appendix A). This decomposition is useful for extracting the corresponding 
scalar functions from the result of a calculation. Although bosonic symmetry of the vertex 
puts some conditions on the corresponding scalar functions, the explicit symmetry of the 
expression is broken, because one of the momenta was substituted in terms of two others. 

To avoid this, one can use a more symmetric decomposition of the general three-gluon 
vertex, proposed by Ball and Chiu fTSUp], 

rMiA.2M3(Pi>P2,P3) = MPi^pI'^pI) ^^^2(^1-^2)^3 +B{pI,pI;pI) ^/.i/.2(pi +P2)m3 
-C{pI,pI;pI) {{PlP2)g^,,f,2 - Pl^,2P^^^^) (Pi - ^2)^3 

+ 1 S(j)IpIp1) {pif,,P2f,,P3f,^ + Plt,2P2t,;P3t,,) 

^Although the standard QCD Lagrangian corresponds to the SU(3) group, our results are vahd for an 
arbitrary semi-simple gauge group. 

^In fact, also completely symmetric colour structures d'^i'^^as niight be considered, but they do not 
appear in the perturbative calculation of QCD three-point vertices at the one-loop level. 

''Another general decomposition of the three-gluon vertex was considered in ref . [E2| . 
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+H{pI,pI,Ps) [-^miM2 {Pl^,■AP^^P3)-P2f,,{PlP3)) + I {pi^,,P2^,,P3^,^-Pl^,^P2^,,P3f,,^ 

+ { cyclic permutations of ipi,fii), fe, 1^2), (Ps, Ats) } • (2.4) 

Here, the A, C and F functions are symmetric in the first two arguments, the H function 
is totally symmetric, the B function is antisymmetric in the first two arguments, while 
the S function is antisymmetric with respect to interchange of any pair of arguments. 
Note that the contribution containing the F and H functions is totally transverse, i.e. it 
gives zero when contracted with any of pi^^, p2^^ or Ps^^- 

Now, before proceeding further, we introduce some notation. For a quantity X (e.g. 
any of the scalar functions contributing to the propagators or the vertices), we shall denote 
the zero-loop-order contribution as and the one-loop-order contribution as X^^\ In 
this paper, as a rule, 

=x(i'^)+X(i'''), (2.5) 

where X*^^'^^ denotes a contribution of gluon and ghost loops in a general covariant gauge 
( |2.3| ) (in particular, corresponds to the Feynman gauge, ^ = 0), while X*^^'^^ repre- 

sents the contribution of the quark loops. 

For example, from (|2.1| ) one can see that at the "zero-loop" level all the scalar functions 
involved in ( p. 4]) vanish, except the A function which is 

= 1. (2.6) 

In what follows, we shall also need to use some other QCD Green functions, including 
those involving the Faddeev-Popov ghosts. As in ( |2.4[ ) we define the corresponding scalar 
structures following the notation of the paper |T^. 

The gluon polarization operator is defined as 

ip) = -S'"'"' {P'9,.,. - P,.P,.) J{p'). (2.7) 
while the ghost self energy is 

The lowest-order results are J'-*^^ = G'-^-* = 1. The one-loop contributions to ^"^^^^{p) and 
n"i'^2(p2) are presented in fig. 2 and can easily be calculated. The results (in arbitrary 
space-time dimension) can be found e.g. in ref. For completeness, we collect the 
relevant formulae in Appendix C. 

The ghost-gluon vertex can be represented as 

f7r'{PuP2;P3) = -ig r"'"' Pi" f,,,ipuP2;P3), (2.9) 

where pi is the out-ghost momentum, p2 is the in-ghost momentum, p^ and fi^ are the 
momentum and the Lorentz index of the gluon (all momenta are ingoing). For F^^g we 
adopt the following decomposition, also used in 



\tJiz{Pl,P2]P3) = 9^l^l^a{P3^P2,Pl) - P3^,P2f,,b{p3,P2,Pl) + Plf,P3^^AP3^P2,Pl) 

+P3^Plt,./{P3,P2,Pl) + Pl^Plf,.AP3^P^'P^)- (2-10) 
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At the "zero-loop" level, 

^^"L=9,,s, (2.11) 

and therefore all the scalar functions involved in ( p.lO| ) vanish at this order, except one, 
a(°) = 1. We shall also need the one- loop-order results for the ghost-gluon vertex in 
arbitrary gauge (the corresponding diagrams are presented in fig. 3). We have calculated 
one-loop contributions to all scalar functions occurring on the r.h.s. of eq. ( p.lOj ), they are 
presented in Appendix D. 

We need F^^g with two Lorentz indices, because this is what enters the Ward-Slavnov- 
Taylor identity for the three-gluon vertex, which, in the covariant gauge, has the following 
form (see e.g. in B, IT^ 



Pf r/,,^2/,3(pi,P2,P3) = -J{pI) G{pI) (g^^l"'' pl-pif,^ pi^'') T^^^^{pi,p^-p2) 

+J{pI)G{pI) (g^rpl-P^^.P^"') ^,,,Ap2,P3;pi)- (2.12) 



It is easy to see that the F and H functions from the three-gluon vertex ( |2.4| ), as well 
as the c and e functions from the ghost-gluon vertex ( 2.1CI| ) do not contribute to this 



identity. Below, we are going to use eq. ( p.l2|) as a non-trivial check on the results for the 



longitudinal part of the three-gluon vertex. 

To conclude this section, we would like to present the notation we use for the integrals 
occurring in the one-loop calculations. We define the integral corresponding to the triangle 
diagram as 

"'^ " / ((p.-g)^)-((p'+g)^)-(g^)-' ^'-''^ 

where n = 4 — 25 is the space-time dimension (in the framework of dimensional regulariza- 
tion|, [^]). A brief overview of relevant results for such integrals in n dimensions is 
presented in Appendix B. It should be noted that all such integrals occurring in the 
present calculation can be algebraically reduced to one non-trivial integral, 

J{l,l,l)=i7c-/'rj^{plpl,pl), (2.14) 

where '^{Pi,P2jPs) = (fi is a. totally symmetric function (see Appendix B for details), and 
three two-point integrals, J(0,1,1), J(1,0,1) and J(1,1,0), which can be expressed in 
terms of a power-like function 

- - - - („-3)(„-4) (--?''"""^ - KT^) ("^) 

as, e.g., 

JjX 1, 0) = i7r"/2 ^ (2.16) 

®For simplicity, we put the dimensional- regularization scale ^ldr — 1- Otherwise, all one-loop ex- 
pressions for dimensionally-regularized quantities should have been multiplied by {/j.dr)^'^- In the final 
results, expanded around n = 4 and renormalized, this scale can easily be restored by inserting ^idb. in all 
non-dimensionless arguments of the logarithms, in order to make them dimensionless. See also section 3.5 
where the renormalization is discussed. 
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and similarly for J(0, 1, 1) and J(1,0,1), where, instead of n{pl) = ^3, we should use 
h{pI) = Ki and k.{pI) = 1^2, respectively. In eqs. ( p.l4| ) and ( |2.16| ), 77 denotes a factor 
constructed of F functions, 

"-T|fr3rr(3-f)^^r(i-H.). (2.17) 
3 Off-shell results 

The set of Feynman diagrams yielding one-loop contributions to the three-gluon vertex is 
presented in Fig. 1. 

When calculating the diagrams, we used the standard technique of tensor decomposi- 
tion |]2^P|, reducing the result to combinations of scalar integrals multiplying the tensor 



structures constructed from the external momenta (see Appendix A). In the Feynman 
gauge, the basic set of scalar integrals (|2.13|) includes the four integrals J(l, 1, 1), J(0, 1, 1), 
J(l, 0, 1) and J(l, 1, 0) only, since massless integrals with two non-positive powers z/j van- 
ish in dimensional regularization [^. For arbitrary ^, we also get integrals with some 
of the powers of the denominators equal to two, see eq. ( |2.3| ). However, with the help of 
the integration-by-parts technique ||2^ these integrals can be algebraically reduced to the 
above basic set (see ref. [^]). While performing the calculations, the REDUCE system 
[ pQ] was heavily employed. 

Before presenting the results, let us define two totally symmetric combinations of the 
invariants formed from the external momenta, 

Q = + {pm) + (P2P3) = -\{pI+pI+pI). (3.1) 



^ = pIpI - {pmf = pIpI - {pmf = pIpI - {.P2P?,Y 

= {PlP2 ) {P1P3 ) + {PlP2) {P2P3 ) + {PlPS ) {P2P3 ) 

= -i MY + {pI? + {vl? - ^pIpI - 2pIpI - 2plpl) . (3.2) 



From the last line of eq. (|3.2|) , one can recognize the structure — 4/C as the Kallen function 
of p\,P2 and P3, see e.g. in fST 



3.1 Results in the Feynman gauge 

Let us consider the one-loop contributions to the three-gluon vertex ( p.4|) in the Feyn- 
man gauge (^ = 0), without the quark loops (the results for the latter are presented in 
Section 3.3). We shall use the standard notation Ca for the Casimir constant, 

jacdfbcd ^ (j^ ^ab ^(j^ ^ ^ ^j^g ^\]{N) group), (3.3) 
whereas the factor t] occurring in the results is defined by eq. ( p.l7| ). 



^An alternative way to decompose triangle integrals (2.13) with tensor numerators was used in pS] 
It was based on a formula from Efll. 
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The one-loop results for the scalar functions 
dimension n, are 



for arbitrary value of the space-time 



1 



b^'^'\pIpIpI) 



(47r)"/2 4(n - 1)/C 

X {{n - 1) (pI + 3(piP2)) [pI{PiP2)V + {Pm)l^l + (P2P3)'«2 +P3'«3 

+4(n - 1)/C [(piP2)v' + ^^?^ - (3n - 2)/C [n^ + /ta]} , 
1 



(3.4) 



^2 



xi (n-1) [(piP3)(P2P3)v5 + (piP3)/«i + {.P2Pz)f^2+pli^^ + (4n-3)/C — ^ 

I Pl~P2 



(3.5) 



1 



X <^3(n-ll 



(47r)"/2 4(n - 1);C 

pI{PiP2)V + (PlP3)'tl + {P2P3)l^2 + pIi^3 - 2 (4:71-3) K. 



Hi — K.2 

Pl-pl 



f^'''\pIpIpI) = v^Ca 



s^'^'\pIpIpI) = o, 

1 



(3.6) 
(3.7) 



X 2 



(47r)"/2 4(n- 1)/C3 
- l)(piP2)(PlP3)(P2P3) + 2(n - 2)p3/C - (n - 7){piP2)}C 

X [P3(P1P2)V5 + (PlP3)'«l + iP2P3)f^2+plf^3 

+2/C [(n + l)(n - 4)(piP3)(p2P3) - (5n - 11)/C] [(piPa)'/' + 
+2p2/C [(n + l)ipiP3)iP2P3) + in- 3)/C] ^ + (4n - 7)/C2 [k, + 



-/C 



2(n + l){piP2){pl - Plf + (4n - 3)/C(p2+p2_2(p,p2))' 

1 



^1 - /^2 I 

Pi - pi J ' 



(3.^ 



h'''\pIpIp^ = 7^.Ca 



(47r)"/2 2(n- 1)/C3 

X {(n^ - l){piP2){PlP3){P2P3) 

X [(PlP2)(PlP3)(P2P3)V5 + (PlP2)(PlP3)/«l + {.PlP2){.P2P3)l^2 + (P1P3) (P2P3) '^3] 

-3(n - l){piP2){PiP3){P2P3)^ [Q^ + Ki + K2 + K3] + 2{n - l)K? if 

+ (n - 2)/C (pi(p2P3) + {PlP2){piP3)) t^l+pl (pliPm) + {P1P2){P2P3)) f^2 

+pI {pliPm) + {PiP3){P2P3)) /ts] } • (3.9) 

When expanded around n = 4, these formulae coincidef^ with the results presented in 
13]. We shall see that the result = is valid also in an arbitrary gauge. It should be 



noted that presenting the results in arbitrary dimension does not spoil their compactness, 
as compared with the formulae expanded around n = 4. 

^"Up to the definition of the renormahzation scheme constant C in which we find to be 

C = — 7 — ln7r + 2 rather than C = —7 — Invr. 
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3.2 Results in arbitrary covariant gauge 



In an arbitrary gauge, the results for the scalar functions of the three-gluon vertex ( p.4|) 
are obviously less compact than those in the Feynman gauge. We list them below, also 
for arbitrary value of the space-time dimension: 



X { [pIpIIC ((8 -A^-{n-2){n- 3)^) pj + 2 (l2 + A{n - 3)^ + {n - 3)^) (pm) 
-ii{{n-A)^+A)ICQ{{n-3){p,p2)Q-{n-A)IC)+a{n-3)^ + 2){n-l)plplpl{p,p2)Q 



X 



-/C 



pI{PiP2)V + {PiP3)k-1 + {P2Pz)t^2+plt^-l. 

((n - 4)e + 4) /C (((n - 4)e - S)p\pI + {{n-2)pl - 2(71- 3) (p.p^))) 

((^-3)^ + 2) in-2)plplplQ] [{piP2)^ + Ks] 
+/C ^ [e ((n - 4)e + 4) /C {{2n - 7)p\plpl + Q {p\{pm) + ^^(^2^3))) 

H {{n - 3)e + 2)pIpIpI (pIQ - 2{n - 4)/C 



/C 



n — 1 



pIpIK (8(3n - 2) + 4(n - l)(5n - 17)^ - 3{n - l){n - A)^' 



H {{n - 4)e + 4) - \)KQ' + e((n - 3)^ + 2)(n - l^p^p^Q [k^ + M ,(3.10) 



{Pl-P"^ 



(47r)"/2 32/C2 p\ pI pI 
{ [-P?P2P3^ (8 - 12^ - + 2)(n - 3)^2) 
+e ((n - 4)e + 4) /CQ ((n - 2)/C + (n - 3)(pip2)b^ + (P1P2))) 
+e ((n - 3)e + 2)p2p2 + (n - l)p^(piP2)(p^ + (^1^2)))^ 

X \^\{p\P2)^ + (PlP3)/«l + (P2P3)/«2 + p1/«3 

-(n - 4)/C [e ((n - 4)e + 4) /C (p^Q + ^2^2^ + ^ ((ri - 3)^ + 2) p?p^p^(p2 ^ (p^^^)) 

X {{PVPi)^ + /«3] 

-/Cp^y, [p2p2^ (8 + 4(n - 5)e - {3n - 10)^') 

+e ((n - 4)e + 4) /CQ2 - e ((n - 3)e + 2) p\plpl{p\ + (^1^2)) 

/C Ki — K2 



1 2 2 

n — 1 P1—P2 



pIpIpIK. (8(4n - 3) + 4(n - l)(5n - 19)^ - [n - l)(5n - 18)^2] 

+e ((n - 4)e + 4) (n - 1)/CQ (p?(piP3) + P2(P2P3)) 

{{n - 3)e + 2) - l)p?p^(p2 _ ^2)2(^2 ^ ^^^p^))] I ^(3^^^) 



c^'^'\pIpIpI) 



^2 r/ 



(47r)"/2 ^ 16/C2 p2 p2 p2 
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< { [^pIpIpIIC (6 + {2n - 5)e + (n - 3)^') 

+e ((^-4)e + 4) QIC (f + {n-Z)pl{pm)) + e ((n-3)e + 2) (n- 1)^2^2(^2 ^ 2 ^^^^^y 

X [pliPm)^ + (PlPs)/^! + {P2Pz)i^2+pIi^Z 

+(n-4)/Cp^ [^((^_4)^+4)/CQ + e((n-3)e+2)p?p^p^] [(piPs)^^ + ^3] 
+/Cp^ ^ [e ((n - 4)e + 4) X: - Q') + e ((n - 3)e + 2)p?p^(p^)^; 

~ [2piP2P3^ (4(4n-3) + 2(n-l)(5n-18)e - (n-l)(2n-7)e') 



n - 1 pf-p^ 



H {{n - m + 4) (n - 1)/CQ (p?(piP3) +^2(^2^3)) 

-e ((n - 3)e + 2) (n - - plf] } , (3.12) 

^^''^^(P?,P2,P3) = 0, (3.13) 



F''''\plphpl)-7^,Cj, 



X 2 



(47r)"/2 32(n - 1)10 pjplpl 
PipIpI [{PiP2){piPz){p2Pi){n - 1) 

X (8(n + 1) + 8(n - 3)^ + (3n^ - 38n + 63)^^ - (n - 3)(7n - 13)^ 
-2/CQ (8(n - 2) + 8(n - 1){ + (n - l)(3n - 14)^^ _ 2{n - l)(n - 3)^^) 
+^{PiP2) (8(n + 3) - (n - l)(n2 - 4n + 23)^^ - 7(n - l)(n - 3)^^)" 
-i {{n - A)i + 4) (n - 1) [((n - 3)^ + 2)(n - ?>){pW2f {pmf {P2Pzf 

+ ((n - 3)^ + n - 1) A:(piP2)^(PiP3)(P2P3) 
- {{n - + 3n - 7) lCp\pl + 3(n - 4);C=^]" 

pI{PiP2)V + (PlP3)«l + (P2P3)«2 +P3'*3 



X 
+2/C 

X [(piP3)(p2P3)(n-4) (8(n+l) + 8(n-l)e + (n-l)(3n-23)e' - 5(n-l)(n-3)e 

-/C (8(5n-ll) +4(n-l)(5n-ll)^ + (n-l)(n^+4n-10)^^ + n(n-l)(n-3)^^) 
-e((^-4)e + 4) (n-l)(n-4) 



X 



{{n - 3)^ + 2) ipiP2)ipiPzf{P2Pzf + A^(PiP2)(piP3)(P2P3) - 3A:^(piP2) 

X [(PlP2)</7 + «3] 

+2/C(/? [p?P2(pi)^ [■'C(n - 1) (16 + 16^ + (5n - 32)^^ - 6(n - 3)^^) 

+^3(^1^2) (8(n + l) + 8(n-l)e+(n-l)(3n-23)C'-5(n-l)(n-3){=') 
-e ((n - 4)e + 4) (n - l)(piP3)(p2P3) 

X an-3)e+2)p2(pip2)(piP3)(p2P3) +((n-3)e + 3)/Cp2(pip2)+^(n-4)/C2^ 

+X:^p?p^p^ (8(4n-7)+4(n-l)(5n-ll)^+(n-l)(13n-30)^2^2(n-l)(n-2)(n-3)^3) 
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2 2 2 
P1P2P3 



X [ki + /t2] 

/C(p^-4(pip2)) (8(4n-3)+4(n-l)(5n-17)^-(n-l)(3n-10)C^) 



X 



+2(piP2)(p? - plYin + 1) (8 - (n - 1)^^ ((^ _ 3)^ + 3) 
2e((n-4)e + 4) (n-l)/C 
IC{pl{pm)+pl{p2P3))-{{n-3)^ + 2){p,p3){p2P3){pl-pl?]] |,(3.14) 



^2 r/ 



C 



(47r)"/2 16/C3p2p2p2 



{[P?P2P3(P1P2)(P1P3)(P2P3) (8(n + l) + 14(n-l)e+4(n-l)(n-7)e'-5(n-l)(n-3)f 
((n - 4)e + 4) Un - 1) ((n - 3)e + 3) (P1P2)'(P1P3)'(P2P3)' + S{n - 4)/C= 



1)^(2 



X [(PlP2)(PlP3)(P2P3)V5 + (PlP2)(PlP3)/«l + {PlP2){P2P3)l^2 + (pm) {P2P3) 1^3] 

-piplpl{PiP2){piP3){P2P3) (24 + 38^ + 4(3n - - 9(n - 3)^=^^ 



+3^ ((^-4)^+4) ((n-3)e+3) {pm)\piP3)'\p2P3f + i{2-i){plf{plf{p''' 



/C 



n 



X [Qip + KI + K2 + «3] 

P?P2P3 (l6 + 4(3n - 8)^ -{n- 2)e - [n - 2)^) 

-e' ((n - 4)e + 4) (n - 4)(pip2)(piP3)(p2P3) 
'^pIpIpI (4(n - 2) + 6(n - 1)^ + {n - l)(2n - 9)^^ - (n - l)(n - 3)^ 
((n - 4)e + 4) ((n - 3)e + 3) (n - l)(piP2)(piP3)(p2P3)] 

(P?(P2P3) + (P1P2)(P1P3)) /«1 +P2 (P2(PlP3) + (P1P2) (^2^3)) /«2 



+P3 (^3(^1^2) + {PlP3){P2P3)) 1^3] } • (3.15) 



One of the main technical problems we met in this calculation was how to bring 
the results for arbitrary ^ to a reasonably short form. Original REDUCE output for the 
numerators of the scalar functions ( p.4| ) was really huge. Then, the problem was how to 
organize the result and which bases to choose. The first basis we needed was one in the 
"space" of the scalar functions if and Ki. It was possible to get better factorization of the 
coefficients by considering not these functions themselves but certain linear combinations. 
Moreover, not all "convenient" combinations happened to be the same for the different 
functions, see eqs. ( p.l0| )-( p.l5|) . Then, the coefficients multiplying these combinations 
of if and Ki are polynomials in C,, n and the momentum invariants. Trying to write the 
latter only as ^1,^25^35 were still getting rather long expressions. The next idea was to 
try to use in some cases also the scalar products {P1P2), {P1P3) and {P2P3), together with 
the notation ( p.l| ) and ( p.2| ) for symmetric combinations. These tricks (as well as looking 
for proper combinations of ^ and n) allowed us to write the expressions in a much shorter 
form. However, this part of the work could not be completely automatized since and 
the scalar products (piPj) are linearly dependent. 
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There are some special values of the gauge parameter ^ we would like to point out. 
First of all, we see that the terms containing Pi,P2 or p| in the denominator disappear 
not only if we put ^ = (Feynman gauge), but also in a "singular" (in four dimensions) 
gauge^^, ^ = — 4/(n — 4). Having no pf in the denominator is especially convenient 
when one considers on-shell limits, i.e., when some of the external momenta squared 
vanish; otherwise, one needs to expand the scalar integrals in the vanishing momenta 
squared (see Section 4). Secondly, many terms vanish for ^ = —2/{n — 3), which could 
be considered an n-dimensional generalization of the Fried- Yennie gauge I^T] (see also 



ref. 1^ and Appendix D). 



3.3 Contributions of the quark loops 

Let us consider the quark loop contributions to the functions ( |2.4| ) . We assume that there 
are Nf quarks which are all massless, and we define 

T^ = lTr(J) = | (if Tr(/) = 4), (3.16) 

where / is the "unity" in the space of Dirac matrices. 

The quark loop contributions do not depend on ^. The results of the calculation are: 

A^'^'KpIpIpD = 7^ N,Tn ^ + , (3.17) 
B^''^\pIpI;pI) = NjTn ^ [k, - n,] , (3.18) 

C^'^'\pI,pI;pI) = 2N,Tn ^ (3.19) 

(47r)"/^ n — I Pi — P2 

S^'^'\pIpIpI)=0, (3.20) 
F^''^\pl,plpl) = --i^N,T^ 1 





V 


(47r) 


n/2 








V 


(4vr) 


n/2 






g' 


7] 


(4vr 


)"/2 



(47r)"/2 ^ " (n- l)(n-2)/C3 
x{2 [{n^ - l){piP2){piPz){P2Pz) + 2(n - 2)pllC - {n - 7){piP2)lC 

X [pliPm)^ + (PlP3)/tl + (P2P3)«:2 +P3'«3 

+2{n + l)(n - 4)/C ^3(^1^2) [{PiP2)v + ^^s] 
+2p2/C ^ [{n + l){pm){p2P^) + {n- 3)/C] + n{n - A)lC [k^ + K2] 

+/C [2{n + l){pm){pl - pI? + {n- 2flC U + pj - 2(^1^2))] | , (3.21) 

P1-P2 ^ V /J I 



H^^>i)ir? T? - ^ ' 2NfTn = 

H [p„P2,Ps) - ^4^^„/2 ^^^f^R _ _ 2)^3 

^^It is not clear whether the second choice could be of use in realistic calculations, since singularity of 
^ in four dimensions requires extra care in renormalizing, etc. 
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X - 1){P1P2){P1P3){P2P3) 

X [iPlP2)iPlP3)iP2P3)^ + {PlP2){Pm)l^l + iPlP2)iP2P3)K-2 + (PlPs) (P2P3) /ts] 

~3{n - l){piP2)iPiP3){P2P3)^ [Q^ + Ki + K2 + K3] + (n - l)(n - 2) /C^ 

+ (n - 2)/C {pI{P2P3) + {PlP2){PlP3)) f^l+Pl (pUPiPs) + {PlP2){P2P3)) f^2 

+P3 (PslPiPa) + {PiP3){P2P3)) i^3\ } • (3.22) 
3.4 Symmetric limit and comparison 



Now, we would like compare our results with those by Celmaster and Gonsalves [|IT] . 
For their study of renormalization-prescription dependence of Green's functions, they 
evaluated the three-gluon vertex function to one loop at the symmetric point, 

pI=pI=pI=p' = -M\ (3.23) 

In this case, we get (P1P2) = (piPs) = (P2P3) = '~\p'^ = ■ At the symmetric point, 
the vertex function simplifies considerably. First of all, because of their antisymmetry, 
the B and S functions ( p.4|) must be zero, 

B{p\p'-p') = S{p\p\p') = Q. (3.24) 

Furthermore, in this limit, the number of independent tensor combinations in the three- 
gluon vertex reduces to three, and the vertex function can be written, in the notation 
used in ]n[, asQ 

rwM2/.3(Pl.P2,P3) = Go{p^) [^^1^2(^1-^2)^3 + ^3(^2 -P3)m + Wi(P3 -Pl)/.2] 

-Gi(/)(p2 -P3)mi(P3 -Pi)m2(Pi -P2)m3 

+G'2(p^) {pif,,P2f,,P3f,^ - Pl^,2P^^,sP3f,,) , (3.25) 



with the three Gi functions related to the scalar functions in ( p. 41) through 

Goip') = A{p\p';p') + \p^C{p\p^-y) + \{p^?F{p\p^-y) + \p^H{p\p\p^),{?,.2Q) 
Gi{p') = C{p\p';p')+'^p^F{p^p^;p'), (3.27) 
G,{p') = C{p\p^-p^) + \p^F{p\p^-y)+H{p\p\p'). (3.28) 

We note that two of these relations may be expressed more compactly as 

G2{p') = G.ip') + H{p\p\p'), Go{p') = Aip'y-y) + \p^G2{p^). (3.29) 

From our results, we obtain the following expressions for the Gi functions (in arbitrary 
gauge and dimension): 

^0 ^^)-J^2^^^^ 

X {pV [s + 12^(14n - 51) + ^i^{n^ - 18n + 60) - (^{n - 4)(n - 12) 



-6/t 



32 + 36^(2n - 7) + 6^^(72 - 4) (n - 6) - i^{n - 4)(n - 3)] } , (3.30) 



^^The counterterm contribution is omitted. 
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{ 



1 



—6ft- 



n 



(47r)"/2 3456 p"^ 
64(n - 20) - 144^(7n - 26) - 24^^(72^ - ISra + 50) - 24n^ + 200ra - 384) 

-4 



1 



64(n-2)-144^(n-l)-12^^(n-l)(2n-7)+r(n-l)(n-3)(n-20) \ , (3.31) 



(47r)"/2 144 p2 



X {pV [128 + 6^(29n-114) + 6^^(^2-1871+60) - ^%n-A){n-12) 



—Qk- 



n - 4 



n 



g^ r] 



8 + 30^(72-1) + 6^2(r2-l)(r2-5) -^^(n-l)(n-3) 

2(3n- 



g'^ V 

(4vr)"/5 

g^ V 

(47r)"/2 



NfT, 



fJ-R 



NfTn 



9{n - 2) 
4 



{2/ ^-3 k] 
4 



4p 9? + 3k 



(4vr)"/2 



9(n - 2) p2 



27 p2 I ' ' n — 1 . 

f 12 — 4 1 
|2p2 (3/7 - 8) - 3fi: A , 



(3.32) 

(3.33) 
(3.34) 
(3.35) 



where (p = f{p'^,p'^,p'^) and k, = K,{p'^). Expanding these results around n = 4 and keeping 
the divergent and finite (in e = (4 — n)/2) terms only, we arrive at exactly the same results 
as Celmaster and Gonsalves, see eqs. (14a), (14b) and (14c) of ||1 l[|p^. 

In ref. ||3^ the QCD renormalization has been considered at an asymmetric point. 



2 2 _ 2 

Pl=P2=P 



pI 



(3.36) 



In particular, the three-gluon vertex was studied at this point, including quark loop 
contributions (with massive quarks). The tensor structures used to decompose the three- 



gluon vertex are presented in eqs. (2) and (4) of ||33[. Of the scalar functions multiplying 



the seven tensor structures defined by eq. (5) of |Q, an explicit result is presented for 
the function Fq only, see eq. (6) of In terms of the scalar functions ( |2.4| ), one finds 

the correspondence Fq ^ —Aij^.Azp^^p^) — B(j? ,Azp^]p^). 

Calculating this combination of the A and B functions for the case ( ^.36| ), we find 
coincidence of the ^ and .^^ contributions[^ as well as the quark loop contibutions in the 
massless limit. However, our result for the Feynman-gauge part is different. So, we do 
not confirm eq. (6) of |^ . 



-'^^Their transcendental constant / is nothing but our p'^'-p{p'^ ^p^ ,p'^)\^_^ (see eq. (B.IO) ) which can be 
expressed in terms of Clausen's function as -^Cl2 (f). 

We note some misprints in [Q: (i) in the definition of the tp function (eq. (8)), x in the denominator 
of the expression under the square root should read z; (ii) before eq. (7), the definition of 1/e should read 
1/e = 2/(Z3 — 4) + 7 + ln(Q^/ {^iiv'^)) (the sign before the logarithm should be changed). 

^^In b is the same as our ^. The integral /(O, z) from is related, in the limit ( |3.36| ), to our 
J(l, 1, 1) as /(O, z) = -{iTr^y^lVPj{l, 1, 1). 



15 



3.5 Renormalization 



In the limit n 4 (e — > 0), the only function which may have an ultraviolet singularity 
is the A function, since this is the only function which does not vanish at the "zero-loop" 
level, see eq. (|2.6| ). In arbitrary gauge, the ultraviolet-singular part of the A function 
follows from eq. ( |3.10|) , 



(4vr) 



-Ca (i + 



R 



K 



(UV) 



(3.37) 



where k*-'^^^ = l/e + . . .is the ultraviolet-singular part of the k function ( |2.15|) . Therefore, 
the divergent part of the counterterm contribution should be equal to minus the r.h.s. of 
eq. ( p. 37 ). This counterterm contribution can be written as[^ 
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(An) 

9^ V 
(47r)2-s 



CaUU 



\NfTn 



+ R 



1 + 



+ 0{e) 



(3.38) 



whereas 
Such a 



where i? is a constant corresponding to the choice of renormalization scheme, 
g"^ = g'^e~"'^{4:Tiy = g"^ exp{e{—'-f + ln(47r))} is the "rescaled" coupling constant, 
re-definition of g"^ is usually performed in the context of the MS renormalization scheme 
which corresponds to the choice R = (for R = 0, eq. ( |3.38| ) corresponds to eq. (15) of 
ref . finil ) . The second line of eq. ( |3.38 ) is more convenient for dealing with the expressions 



obtained in the present paper, since one can keep g'^rf/^ATi)^^'^ as an overall factor. Here, 
we have used the fact that 



-7£ 



(1 



12" ^ 



+ o(r 



(3.39) 



The term in ( p.39 ) is not relevant for the ultraviolet renormalization at the one-loop 
level. However, it yields finite contributions when one has l/e^ infrared (on- shell) singu- 
larities (see section 4.3). 

If we now recall the existence of the dimensional- regularization scale parameter 
(which we usually put equal to one, see footnote at the end of section 2), we shall see 
that, as opposed to other one- loop-order contributions, the counterterm (|3.38|) should not 
be multiplied by the factor (/iDi?)^^, 
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see also in ref. 
Examining eq. (|3.37|) , it is interesting to note that for 



71=4 



Ca 



(3.40) 



we do not have any ultraviolet singularity in the one-loop contribution to the three-gluon 
vertex. If we substitute Nj = 6, Ca = 3 and Tr = i, we get 



4 On-shell limits 

There are two main on-shell cases of interest: when one or two of the external momenta 
squared are zero. However, it is also instructive to consider, as a separate case, the limit 
when one external momentum (not only its square) vanishes. 

^^In commonly-used notation, the expression ( |3.38| ) corresponds to {Zi — 1) at the one-loop order. 
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4.1 One external momentum squared is zero, pl = 

In this case, we should substitute in the expressions for all scalar functions 



0, {P1P2) = -UpI + pD, (pm) = -UpI-pI), (P2P3) = -^2)- 



Note that now we should consider the scalar functions A,B,C and F from 
permuted arguments as well. 

The result for the triangle integral ( p.l4|) simplifies in this limit, 



(4.1) 
with 



J(l,l,l) 



ITT 



2-e 



V f{pi,pio) 



P3-- 



-pI)-' 



-pI)-' 



pi 



pI 



(4.2) 



where 77 is defined by eq. ( p. 17] ). Moreover, in the framework of dimensional regularization, 

J(1,1,0) 



23 



0, 



(4.3) 



while the results for J(1,0, 1) and J(0, 1, 1) remain unchanged. As to the factor l/e^ in 
( ^.2| ), one power of e is cancelled by the expansion of the numerator, while another power 
of e survives and corresponds to the infrared (on-shell) singularity which arises in the 
scalar integral in the limit p| = 0. 

For the Feynman gauge, = ^ (and also for the singular gauge, ^ = — 4/(n — 4)), it is 
enough to perform the above substitutions to get the answer. In the case of arbitrary ^, 
however, the situation is more tricky, due to the presence of p\ in the denominators of the 
scalar functions. Here, in order to get a correct answer, one needs to expand the integral 
J(l, 1, 1) in P3 and keep the term of order p|. 



■/(1, 1,1; 



n- 



i tt" ' ^ - 



-pIY' 



-pI)-' 



Pi 



pI 



:i+e){pi-pi) 



2\2 



^-e)[{-plr + {-pi) 



pI 



-Plf-'' 



-pIY-' 



pI 



pI 



-0({plf).{AA) 



To present the results obtained for the scalar functions 
nient to introduce 



5 



pI 



pI 



12 



Pi 



pI 



in this limit, it is conve- 



(4.5) 



In this section, we present such results for the three-gluon scalar functions in the Feynman 
gauge, and also for the quark loop contributions. The expressions for arbitrary covariant 
gauge are listed in Appendix E. We present the results for the A,B,C and F functions 
of the arguments (Pi,P2;0) i^^Phpl)- The results for the third set of arguments 
corresponding to the cyclic permutations in eq. ( |2.4| ), {p2,0]Pi), can be obtained from 
the functions of {Ojpf^p'^) by using the symmetry (for the A, C and F functions) or anti- 
symmetry (for the B function) with respect to the first two arguments, and interchanging 
Pi ^ P2- The H function with permuted arguments does not change, while the S function 
is zero (at one loop). 
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The results for the gluon and ghost contributions to the three-gluon scalar functions 
in the Feynman gauge are: 



(47r)"/2 4(n- l)(n-4) 

X {(n - 4)(3n - 2) [k^ + K2] - n{n - l)(5i2)~^ [«i - ^2]} , (4.6) 



(47r)"/2 4(n - l)(n - 4) (pf - p^) 
X {{n - 4)fi;i ((2n - l)pl - {An - ?,)pt) - 2(n - l)(5i2)-' [«! - «;2] (2p? - 3p^)} , (4.7) 



(47r)"/2 4(n- l)(n-4) 



(4/1^ - 21n + 14) [ki - K2] , (4.8) 



(47r)"/2 4(n- l)(n-4) (p2-p2)2 
X {-(n - 4)(p2 -p2)^^ ^^g^ _ 5)p2 _ _ 3^^2 

+2(n - l)p? [k, - K2] {{n - 3)p? - (n - 5)p^) } , 



c^'^'\pIpI^) 



(47r)'*/2 2(n-l) (p?-pi) 



(4.9) 
(4.10) 



(47r)"/2 ^ 2(n - l)(n - 4) pf (pf - pi)2 
x{(n-4)(p? -p^)«i (np? - (4n-3)p^) + 6(n-l)p?p2 [«i - k^]] , (4.11) 



F(^'°)(p?,p^;0) = 



c 



(47r)"/2 ^ (n - l)(n - 4) (p? - pi)2 

x(4n - 7) {(n - 4) [ki + ^2] - 2(5i2)-' K - «2]} , (4.12) 



f^'''\^.pIpI)--^^Ca 



(47r)"/2 " (n - 1) (n - 4) p? (pf - p2)2 
x{(n - 4)(p2 - pI)ki [(3n + 1) - (n + 2)(<5i2)-' - 6(5i2)- 

+2p2 [/ti - K2] [{n^ - I2n + 17) + 4(n - 4)(5i2)"' + 6(^12)-'] } , (4.13) 



h^'''\pIpI^) = -AZ,c. 



(47r)"/2 ^ (n - l)(n - 4) (p? - p2) 
X {3(n - 4)((5i2)-' [ki + K2\ + K - K2] ((n^ - 2n - 2) - 6((5i2)-')} . (4.14) 
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The quark loop contributions in the hmit p| = are: 



9 V 
(47r)"/2 



n-2 



Tr [ki + K2\ 

n — 1 
n-2 



n — 1 
n-2 



c^'''\pIpI;0) 



g'^ 1] n — 2 



9 V 
(4vr; 



n — 1 

irp^ "^NfTfi — ^ 1 



C^''''\0,pl-pl) = j^2NfTn 



n — 2 Ki — K2 
- 1 pI-: 

n — 2 Ki 



(47r)"/2 ^'y^^ n-1 pV 



F^'-\pIpI;0) 



g"^ T] 



ANfTn 



n 



(47r)"/2 (n-l)(n-2) [pl-plY 

X {(n - 4) [ki + K2] - 2((5i2)"^ K - K2]] 



(4.15) 

(4.16) 
(4.17) 
(4.18) 
(4.19) 
(4.20) 



(4.21) 



F(^''')(0,p?;p^) 



^2 r] 



^NfTn 



(47r)"/2 " (n - 1) {n - 2) (n - 4) (p? _ p^) 

x{{n-A)Ki [{n^ - 5n + 8) - {612)'^ {n + 2)- 6(612) 



-2 



+16 



pI pI 



{Pl - plf 



[ki - K2] {{n - 1)pI - (n - A)pt) | , (4.22) 



(47r)"/2 ^ « _ i)(^ _ 2)(n - 4) {pl - pl) 
X {3(n - 4)(5i2)"' [ki + K2] + [(ri - 2)(n=' - 5n + 7) - 6(5i2)"'] [ki - K2]} . (4.23) 

We note that there is an interesting relation between the one-loop contributions to the 
A and B functions of permuted argumentsf^, 

{pI + pI)a^'\pIpI^) - p\a^'\q,pIpI)-pIa^'\^.pIp\) 

+ pI B('\0,pI;pI)+pI B('\0,pI;pI) = 0. (4.24) 

This relation is also satisfied by the expressions for arbitrary ^ given in Appendix E. 

The infrared 1/e singularities of the results for gluon and ghost contributions (also 
in arbitrary gauge, see Appendix E) have been compared with the results given in ||14|| , 
eqs. (24)-(25). The functions defined in are proportional to our functions Zijk (see 



^^This is also valid for the zero-loop functions, since A^^^ = 1 and i?^"^ = 0. 
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Appendix A) which can be represented as hnear combinations of the scalar functions (|2.4|) , 
including those with permuted arguments. To get renormalized results, the counterterm 
(|3.38|) was added to all A functions. In the MS scheme, the obtained results coincide^ 
with those presented in ref. [H, eq. (25). 



4.2 One external momentum is zero, j^s = 

In this case, pi = —p2 = p {Pi = P2 = P^)? and the proper limit of eq. ( [4.2| ) yields: 

1 



J(l,l,l) 



1 TC 



2~e 



P3=0 



2\-l-e 



(4.25) 



Actually, we get some powers of {pf —pi) in the denominator from the /C's, eq. (|3.2|), since 



/C 



^(Pi ~ P2Y in this limit. Therefore, we should be careful taking the limit p^ — > pf 



and expand the numerator up to higher powers of 



5[ 



12 



Pi 



26 



12 



l + S 



(4.26) 
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Starting from the result for the case = 0, we need to express K2 as ki times an expansion 
in 5^2- This can easily be done using 



j(i, 0, 1) = j(o, 1, 1) (4) ' = j(o, 1, 1) f: ^^{6[,y. 



(4.27) 



In practice, we need the terms up to only. It is interesting that in this special case 

we do not need the p"^ term of the expansion of J(l, 1, 1) (like in eq. (|4.4| ) ) since it cancels 
in all contributions. 

In this limit, there are only three independent tensor structures left, and the coefficients 
multiplying these structures can be expressed in terms of the "surviving" scalar functions 
as0 



{p, -p, 0) = 2 g^,,,p,, A{p\p'- 0) +p' C{p\p'- 0) 



- (^mM3PM2 + ^M2/.3Pw) [^(0,p';p') - 5(0,/;/)] - 2p^,p^,p^^ C{p^,p^; 0) . (4.28) 

The one-loop contributions to the scalar functions appearing in ([4.28|) are: 
A(^'«)(/,/.o)+p2 c(i'«(/,/.Qy 



9 V 



(471)^^/2 32(n- 



— {8(3n-4) + 12^(n-l)(n2-4n+ 2) -^^(n-l)(n-4)(n+2)} , (4.29) 



(0, /; /) - (0, /; /) 1 = A^W (p2^ p2. 



9 V 



(47r)"/2 8(n-l) 



^ [2{n^ - 7n + 4) - 2^(n - l)(4n - 13) + - l)(n - 4)} , (4.30) 
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Up to a misprint in eq. (25g) of where, in the term proportional to A (their A is nothing but our 



£,), the contribution + 2p) should read ^(1 + + 2/ 
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This corresponds to the decomposition used in ref. 1 14 , eq. (20) 
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g'^ 7] (n - 4) K 



n 



(47r)"/2 32(n - l)p2 
_ 5(1,9) (0 



2 2-1 



(4vr)"/2 



4^(n-l)(3n-8) - ^2(n-l)(n-2)} ; 

(4.31) 
(4.32) 



5(2 7] 



NfT; 



(n 



fJ^R 



n — 1 



9^ V 
(47r)"/2 



NfT, 



fJ-R 



(47r)"/2 

{n-A){n-2 
n — 1 



^'"^ 2iV;r«^«:, (4.33) 



n 



P 



(4.34) 



where k = K,{p^). 

We note that, according to eqs. ( [4. 301 ) and ( [4.33| ), the following relation holds for the 
zero-momentum case: 



(p2, p2. 0) _ ^(1) (0, p2. p2) ^ ^(1) (Q^ p2. p2 



0. 



(4.35) 



This also follows from eq. ( ^.24| ) . The relation ( |4.35| ) is valid for arbitrary values of n and 
^. Using ([4.35|) , we can reduce the number of tensor structures in ( [4. 281) from three to 
two, namelyp^ 



P 9f^lfi2 



P.^pJ C^'\p',p';0), 



(4.36) 



where the results for the scalar functions are given in eqs. ( [4.30| ), ( [4. 311 ), ( |4.33| ) and ( |4.34] ). 
Note that the first tensor structure on the r.h.s. of ( 4.36| ) coincides with the "zero-loop" 
vertex structure, given by (^.Ij). 

To make the complete comparison of our expressions (|4.29| )-( ^.31| ) (contributing to the 
three-gluon vertex at = 0, eq. ( [4.28|) ) with eq. (20) of [|14|, we need to renormalize our 
coefficients ( [4.29| ) and ( |4.30| ) by adding the counterterm ( p.38| ) to all A functions involved 
and putting n = 4. Performing this in the MS scheme (i.e. at i? = 0), we find that our 
results give the same as eq. (20) of |1^ . 

We have also compared the renormalized (in the MS scheme) version of eq. ( [4.36|) with 
the one-loop results presented in ref. [^, eq. (AlO). According to eq. ( |4.36| ) there should 
be the following correspondence between the functions Ti and T2 used in [^] and A^^^ 
and C(^): 

Ti(/) ^ A«(p^/.o)^ T^ip') ^ -2p'C^'\p',p';0). (4.37) 

Renormalizing A^^\p'^ , p'^ ; 0) (given by the sum of eqs. ( |4.30| ) and ( [4.33|) ) and putting 
n = 4, we arrive at the same result for Ti as given in eq. (AlO) of ref. ^5|]0. The result 
for T2 is finite as — > 4 and should correspond to the sum of our eqs. ( |4.31| ) and ( [4.34|) 

^"This corresponds to the decomposition used in ref. jss], eq. (A2). 

^^The gauge parameter used in |^ corresponds to our 1 — ^. To avoid confusion, we shall call their 
parameter S.bl = Their constant T corresponds to our N jTp. The results presented in p5| are taken 

at = — /i^, where corresponds to our /^Ij^j (see the discussion in section 2). Thus, putting — 
formally corresponds to omitting the terms containing ln(— p^) (which appear due to the expansion of 
k{p^) in e) in the renormalized expressions. 
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taken at n 
[00 



4. However, our result for T2 is different from the one given in eq. (AlO) of 



4.3 Two external momenta squared are zero, Pi = pi = ^ 

In this case, we substitute 

pI=pI = 0, pI = p^ (P1P2) = |p^ (piPa) = {P2P3) = -\p^] 



J(l,l,l) 



i7r2-^r7(^(0,0,p^) 



p2=p2=0 

J(1,0,1) 



■ 2^e 1 



0, 



pf=pi=o 



(4.38) 
(4.39) 
(4.40) 



= ^(0,1,1) 

while the result for J(1,1,0) remains unchanged, eq. ( p.l6|) . Note that now, when two 
external lines are on shell, the infrared singularity in J(l, 1, 1) is stronger and gives 1/e^, 
eq. O). 

Again, it is enough to make the above substitutions to get the result in the Feynman 
gauge (^ = 0), and in the singular gauge (^ = — 4/(n — 4)), but the situation is more 
tricky for arbitrary ^ since we have and P2 in the denominators of the scalar functions. 
To solve this problem, we need to consider the expansion of J (1, 1, 1) in and 

Two independent ways were used to get the results for the scalar functions in this 
limit: 

(i) We take the expressions for one of the momenta squared equal to zero (see Section 4.1), 
and put the second momentum squared equal to zero. In the corresponding expressions 
(see Appendix E), all pj occurring in the denominators are always accompanied by the 
corresponding = i^ipf) in the numerator, which should be put equal to zero when pf 
vanishes. 

First, we put J(l, 0, 1) = J(0, 1, 1) = 0. Since the resulting expressions did not have 



n 



singularities for pi = p\, the next step was to put pi = p^ = p^. Then, the integral 
J(l, 1, 1) was expanded in p\lp^ keeping the terms up to {pl/p^Y ■ And finally, the limit 
was taken. The following formula was used to expand J(l, 1, 1): 



pI 



J(l,l,l) 



1 TT 



2-e 



p\=pI=pI 



-P 



2Fl 



+ e 



' 2 

1 + 26 



p2 



2F1 



+ e 



' 2 
1 + 2e 



E 



;i + 2e), 




+ 



(4.41) 



(4.42) 



The results obtained in these two ways coincide, and the expressions obtained for the 
scalar functions ( ^^ are presented below. Due to the symmetry properties, the functions 
A,C and F of the arguments (0,p^,0) are equal to the corresponding functions of the 

^^Their result for T2 is proportional to [{—^+S,bl)Ca + ^T], whereas our expressions yield 



[(— II + |CsL + dA + So, the quark contribution is the same while the sum of gluon and 

ghost contributions to T2 is different. Note, that in the Feynman gauge {^bl = 1) our results coincide 
and yield [— |Ca + |t] . Thus, the disagreement does not influence the one-loop part of the two-loop 
result presented in the Feynman gauge, eq. (B4) of [Esl. 
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arguments (p^,0,0), while the B function with these arguments permuted changes sign. 
The H function is the same for all permutations. 

The resulting one-loop contributions (without quark loops) to the scalar functions 
(|2.4|) in arbitrary gauge are: 



^^"'^(0' ^-y) = J0r2 Ca ^^^^ {48 - S^n -3){n-6)+ e{n - Af} , (4.43) 

A(i'«)(p2, 0; 0) = --^^ Ca \-, r 

' ' ^ (47r)"/2 ^ 64(n- l)(n-4) 

X {16(2^2 - 13n + 8) + 4^(n - l){2n - 9)(5n - 16) - 5^^(n - l)(n - 4)^} , (4.44) 

5(i'«)(0,0;j92) = 0, (4.45) 



5(^'«(/,0;0) = -/^a 



(47r)"/2 64(n- l)(n-4) 
X {l6(4n2 - 21n + 14) + 4^(n - l)(10n2 - 79n + 152) - 5^2(n - l)(n - 4)^} , (4.46) 

C(^'«)(0,0;/) = ^|-^ 4(n - 4)(n - 6)^^ 



X 



{l2(n - 6) + 2i{n'^ - lln + 36) - i^{n - 3)(n - 8)} , (4.47) 



C(i'«)(p^O;0) = -7^^ — ^7 K 

' ' ' (47r)"/2 ^ 32(r2- l)(n-4)p2 



X 



{l6n{n - 4) + 4^(r2 - l)(6n2 - 41n + 72) - 3^'^{n -l){n- 4)^} , (4.48) 



^ ' '^^ (4^)„/2 2(n-l)(n-4)(n-6)(j92)2 ^ 



X 



{4(n - 3){n - 6)(n - 13) - 2^{n - l){3n^ - 20n + 18) 

+^2(n - l){n^ - 16^2 + 74n - 78) - - l)(n - 3)(n - 8)} , (4.49) 



.2 



F(l'«)(p^O;0) = -,f^a 



(47r)"/2 16(n - l){n - 4)(n - 6)(p2)2 
X {l6(n - 6)^(4n - 7) + 4^(n - l)(n - 6){10n^ - 87n + 152) 

+^^{n - l)(5n3 - 94^2 + 504n - 624) - 8^\n - l)(n - 3)(r2 - 8)} , (4.50) 

H^^'^\0,0,p^) = -j^^ Ca iQ^n - l){n - A){n - Q)p^ 

X |l6(n - 4)(n - 6)(n + 5) + 24^(n - l){n - Q){n^ - 6n + 12) 

-i'^in - l){7n^ - llOn^ + 532n - 696) + 4^^(n - l){n - 3)(n - 8)} . (4.51) 
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The quark loop contributions yield: 



(0,0;p2) = 5(i'^)(0,0;p2) = C(i'«)(0,0;/) = 0, 
(p^ 0; 0) = (p^ 0; 0) = NjTn ^ " ^ 



(P',0;0) 
(0,0; p2) 

(P',0;0) 
(0,0, 









n/2 








V 


(47r) 


n/2 








V 


(4vr) 


n/2 








V 


(47r 


)"/2 



2NfTR 



(47r)"/2 
n — 2 K, 
n — 1 p2 ' 
1 



n — 1 



(r2-l)(n-2) (p2)2^ 
n{n — 6) K 



(n- l)(n-2) (p2)2' 
ri^ — 3?2 + 8 K 



(n — l)(n — 2) p2 



(4.52) 
(4.53) 

(4.54) 

(4.55) 

(4.56) 

(4.57) 



Since in the limit p\ = pi = ^ the scalar functions may depend on p\ = p^ only, the 
independent tensor structures (expressed in terms of pi and P2) can be chosen antisym- 
metric with respect to the permutation (pi, /ii) ^ (p2, Ai2)- The three-gluon vertex can in 
this limit be written as 



^^l,^l2^li{Vl,V2,Pz) 



+ 



^WM2(Pi -P2)m3 Ui{p^) + 



9fJ.lfJ.iPlfA.2 9fJ.2fi3P2fj.i 



U2{P') 



U-iip'^) + Pl,,,P2nAPl - P2)f.a Ui{p^ 



+Pl^2P2^l(pi -^2)^3 ^5(P^) + 



f/6(p' 



+ 



Pl^.^Pl^^2P2^.3 - P2^.,P2^,2P^^,, 



Urip'). 



(4.58) 



This decomposition is analogous to eq. (29) of Ijl^, and the functions Ui are proportional 
to the functions Fj used in |]14| (we have different numbering, see in Appendix F). 

The following representations of the Ui functions in terms of scalar functions corre- 
sponding to the decomposition ( p. 41) can be derived in this limit: 



U,ip') = AiO,0;p') - \ p' 0(0,0; p') - i (p^ F{0,0;p') + \ p' H{0,0,p'), (4.59) 



U2{p' 



-2A(p2,0;0)-/C(p2,0;0), 



(4.60) 



Usip') = -Aip\ 0; 0) -B{p', 0; 0) - i p' C{p\ 0;0) + \ {p^f F{p\ 0; 0) + | p^ H{0, 0, p' 



U,{p') = C{p\0;0)-lp' F(p2,0;0) 
U,ip') = 2 C(p^ 0; 0) + CiO, 0; p') + \ p' F(0, 0; p') - /f(0, 0, p' 



(4.61) 
(4.62) 



(4.63) 

[/6(p2) = 2C(p^0;0), (4.64) 
Urip^) = -C(p^0;0) + ip2^(p^0;0) + i^(0,0,p2). (4.65) 

The explicit results for the Ui functions are presented in Appendix F. The infrared- 
divergent contributions were successfully compared with eq. (30) of [0, where the corre- 
sponding contributions to Fj are presented. Note that f/4, f/5, Uq and f/7 can be directly 
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compared with [Tj], while the expressions for Ui, U2 and U3 (containing the A-function) 
should be renormalized by adding the conterterm contribution ( p.38| ) to all A functions 
involved. 

The result for the three-gluon vertex in the Feynman gauge (for pf = P2 = 0) is 
available in Appendix B of ref. |[15|. It is expanded around n = 4, and the divergent and 
finite (in e) parts are presented. In this limit, our expressions yield the same as the results 
of 0. 



5 Conclusions 

In this paper, we have presented results for the one-loop three-gluon vertex valid for 
arbitrary values of the space-time dimension, n, and the covariant-gauge parameter, ^. We 
have considered the general off-shell case (arbitrary ^1,^25^35 see section 3), as well as all 
on-shell cases of interest (section 4). Moreover, having the results in arbitrary dimension, 
it was possible to get all on-shell expressions just by considering the corresponding limits 
of the general (off-shell) results. This would be impossible if one started from the off-shell 
results expanded around n = 4, because in this case the infrared (on-shell) divergencies 
would appear as logarithms of vanishing momenta squared. The only restriction we used 
in our calculations was that in the quark loop contributions the quarks were taken to be 
massless0. 

To calculate the vertex, we used the decomposition ( |2.4| ) (adopted from the paper |T^) 
and considered the six scalar functions. A, B,C, S, F and H, which completely define the 
three-gluon vertex. One of these functions, namely the S function, was found to be 
identically zero at the one-loop ordeiQ, see eqs. ( |3.13| ) and (|3.2CI|) . We have also checked 



that S = when massive quarks are considered. It is not clear whether it vanishes also 
at the two- loop level. For the five remaining functions, A,B,C,F and H, the general off- 
shell results are given in eqs. ( ^^ -( ^^ (Feynman gauge), eqs. (t3.10|) - (t3.15|) (arbitrary 
gauge) and eqs. ( p.l7] )-( p.22| ) (quark loop contributions). They involve only one non- 



trivial function, '^{Pi,P2,P3) (see eq. ( p.l4|) and Appendix B), which is related to the 
scalar one-loop triangle diagram. For special cases, we have successfully compared our 
results with those from the papers [0, |T^ (for details, see section 3). 

Starting from general expressions and putting some external momenta squared equal 
to zero, we considered the on-shell cases p^ = (section 4.1), ps = (section 4.2) and 
Pi = P2 — ^ (section 4.3). For all these cases, the results in arbitrary gauge and dimension 
were presented, see section 4 and Appendices E and F. For special cases, our results have 



been compared with those presented in refs. [|14| , ^ . Thus, we can see that Table 1 from 



the Introduction is completely filled in! Moreover, all results are valid for an arbitrary 
value of the space-time dimension. Thus, the only thing which at the one-loop level is left 
for the future is to allow for non-zero quark masses in the quark loops. 



Furthermore, we have obtained general results for the ghost-gluon vertex ( p.lO|) , see 
eqs. ( p.7| )-( [D.ll| ) in Appendix D. Employing these results, together with two-point con- 
tributions listed in Appendix C, we have checked that the Ward-Slavnov-Taylor identity 

^^This restriction does not affect any of the results for the gluon and ghost loop contributions (the 
functions marked (1,C)) which are indeed the most general ones. 

In the Feynman gauge and four dimensions, this results was obtained in U3. 
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(|2.12|) for the three-gluon vertex is exactly (i.e. for arbitrary n and ^) satisfied by the 
expressions obtained, as it has to! This was another non-trivial check on the longitudinal 
part of the vertex {A, B, C and S functions). 

We note that techniques have recently become available to study the off-shell 

massless vertices at the two-loop level, at least in the Feynman gauge. Here, the main 
difficulty are the integrals with higher powers of irreducible numerators P7[] . 

In the future, the one-loop quark-gluon vertex can also be considered in a similar 
way. For the case of massless quarks, one can use the same approach as in this paper. For 
massive quarks, one should study in more detail what the corresponding scalar integrals in 
arbitrary dimension are^. We note that some results for the quark-gluon vertex (and also 
for the QED vertex which formally corresponds to one of the two diagrams contributing 



to the quark-gluon vertex) can be found e.g. in |35 



Acknowledgements. The authors are indebted to J.S. Ball and T.-W. Chiu for useful 



communications and, in particular, for confirming two minor misprints in |T3| mentioned 
in section 3.1 and Appendix D. O.T. is grateful to NORDITA for supporting his visit 
to Bergen where this work was initiated, and to the University of Bergen for hospitality. 
A.D.'s and P.O.'s research was supported by the Research Council of Norway. 



25 



This is also the reason why we did not consider massive quark loops in this paper. 
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Appendix A: Decomposition of three-gluon vertex 

If we express in terms of the two other momenta, = —pi — P2, we get the following 



decompostion of the three-gluon vertex ( |2.2| ) 



+fi'/ii/i2P2^3^002 + fl'Mi/i3P2^2^020 + fl'M2M3P2^i^200 
+Pl^,^Pl^,2P^^,3^1U + P2^,,P2f,^P2f,,Z222 

+P1miP1M2:P2;.3^112 +Pl^iP2;.2Pl^3^121 +P2;.iPl^2Pl^3^211 
+Pl^^,P2^,^P2f,.,Zl22 + P2^,Plf,^P2f,3Z212 + P2f,,P2f,^Pl /i3 221 

where Zjki are scalar functions depending on pfyp"^ and p|. 



Comparison with the decomposition ( |2.4| ) used in |]T3| gives the following representa- 



tions of Z's in terms of the functions (|2.4D used by Ball and Chiu 



Zooi = A{pi,pl;pl)-ip,p2)Cipl,pl;pl)+B{pl,pl;pl)HPiP2)iP2P3M 

(A.2) 

^002 = -A{pI,pI; pI)+ {pip2 )C{pI,pI] pI )+ B {pi ,pl;pl )- (pipa ) {piP3 ) ^ (p? , P2 ; Pa )+ (PiPs ) H, 

(A.3) 

(A.4) 

Z200 = 2A{pI,pI;pI) - 2(p2P3)C^(P2,P3;P?) -Pi(P2P3)i^(P2,P3;pD (A.5) 
Zolo = -'2Aiplpl;pl) + 2{p,ps)C{plplpl)+plipiP3)F{plpl,pl)-plH, (A.6) 

^020 = -A{pl,pl;pl'yr {piPs )C{pI,pI, pI )- B {pi ,p\]pl )- {piP2) {pm )F{pI,pI; pj )+ (pips ) 

(A.7) 

Zni = 2Ciplpl,pl)+plFiplpl,pl), (A.8) 

^222 = -2C(p2,p2.p2) _^2^(^2^^2.^2)^ 

^112 = -C(p2,p2.p2) _^ (p,p,)Fiplpl;pl) + H-S, (A.IO) 

^121 = C{pIpI;pI) - {piP2)F{pIpI;pI), (A.ll) 

Z122 = -C{pIpI;pI) + {piP2)F{pIpI;pI), (A.12) 

Z2u = Cipl,pl;pl) + 2Ciplpl;pl)-ip2P3)Fipl,pl;pl)+plFiplpl,pl)-H~S, (A.13) 
^212 = -C'(p?,p^;pD-2C(p^,p2.p2)^^p^p^)^(p2^p2.p2)_p2^(^ (A.14) 

^221 = _ (p,p,)F{pIpI;pI) -H~S, (A. 15) 

where = H{plplpl) and 5 = S{pj,plpl). 

Solving these equations we get the following results for the scalar functions (|2.4|) , 
including those with permuted arguments, in terms of Z's: 

Siplplpl) = \ {-Z112 + ^121 + ^122 - ^221} , (A. 16) 

H{pIpIpI) = \ {Z112 + ^121 - ^122 - ^221} , (A. 17) 
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a{pI,pI,vI) = i I (P1P2) [-^111 + ^222 + ^211 - ^212] + Zqqi - Z002 - {vl + pI)h^ , 

(A.18) 

A{pIpI;pI) = \ {-(P2P3)^222 + ^200 " pIh] , (A.19) 
MpIpIpI) = \ {biP3)^iii + Z^^^-vlE] . (A.20) 

B{pI,pI;pI) = \ \p\ [Z112 - Z122] +pI [-Z121 + Z221] 

~{PlP2) [Zui + ^222 — ^211 — ^212] + ^001 + ^002 + pIs}, (A.21) 

B{pIpI;pI) = I {-(P2P3)(^222 - 2Z122) + ^200 - 2Z,oo + {pI " pI)h} , (A.22) 
B{pIpI;pI) = '2 {-{PiP3){Zin - 2Z121) + ^010 - 2Z020 + {pI-pI)h} , (A.23) 



1 

Pl~P2 



c{pl,pi,pl) = - — 2 {(P1P3) [^111 - ^211 - ^121 + ^221] 



+ (^2^3) [^222 + Zu2 - Zx22 - Z212]} , (A.24) 

C{plpi,pl) = {{PlP2)Z222 +plZl22] , (A.25) 

P2 ~ Ps 

C{pIpI;pI) = {{piP2)Zni +pIZ,2i} , (A.26) 



F{pI,pI;pI) — — 2 {^111 + ^222 + ^112 — ^121 — ^211 " ^122 — ^212 + ^22l} , (A. 27) 

Pi ~ P2 

F[pIpIpI) = ^2^2 {Z222 - 2Z122} , (A.28) 

P2-pi 

HpIpIpD = {Ziii - 2^121} • (A.29) 

p-i ~ Pi 



Appendix B: Scalar integrals 

As mentioned in section 2, the results for the scalar functions occurring in eqs. ( p.4|) and 
( |2.1CI| ) can be represented in terms of the following Feynman integrals, corresponding to 
a scalar one-loop triangle diagram: 

"'^ " / {{P2-qf)H{pl+<iYY^{<i'Y^' ^""-'^ 

where n = 4 — 2^ is the space-time dimension. 

When we perform calculations in the Feynman gauge and express the scalar numerators 
in terms of the denominators, the powers of the denominators, z/j, can be one, zero, or even 
negative. In an arbitrary gauge, the integrals may have powers of Ui equal to two, due to 



the presence of (p ) in the ^ term of the gluon propagator, eq. (|2.3|) . Nevertheless, by 



using the integration-by-parts technique |^ these integrals can be reduced to those with 
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z/'s equal to one or zero. The corresponding algorithm for the integrals ( [B.l| ) is described 
in detail in [00- 

Then, if two or three z/'s are non-positive integers, the dimensionally-regularized inte- 
gral ( [B.l| ) vanishes since it corresponds to a massless tadpole. When one of the z/'s is 
zero, the integral ( |B.1|) corresponds to a two-point function. Therefore, it is proportional 
to a power of the external momentum squared times some F functions with arguments 
involving n and z/'s. We mainly need the result for two remaining z/'s equal to one (e.g. 
J(l, 1, 0)) which is given by eqs. ( |2.15D -( p.l6| ), in which case this F factor is nothing but 
r^, eq. 

Then, the integrals with one negative z/ can easily be reduced to integrals with the 
corresponding v equal to zero, for example: 



J(l,l,-1) 
J(l,l,-2) 

J(l,l,-3) 



-(P1P2) J(i,i,o) 



n — 1 
1 

n — 



2 2 
P1P2 



J(1,1,0), 

n + 2){p^p^f - ?>py^ ip,p2) J(l, 1, 0). 



(B.2) 
(B.3) 

(B.4) 



Thus, the only non-trivial function which occurs in our calculations is (p which is 
related to the triangle integral J(l, 1, 1) (cf. ( |2.14| )) via 



J{1,1,1)= in''/' vM,pIpI), 



(B.5) 



where rj is defined by eq. ( |2.17|) . 

In fact, the general results (i.e. for arbitrary n, z/j and pj) for the integrals (|B.1|) are 
available [|39|. They can be represented in terms of Appell's hypergeometric function of 
two variables, F4. As dimensionless variables, one can use 



X 



and 

pi 



_ pI 
y = — ■ 

pi 



(B.6) 



When z/i = z/2 = z/3 = 1, all the F4 functions can be reduced to 2-^1 Gauss hypergeometric 
functions of more complicated arguments, by using reduction formulae for F4 functions 
(see e.g. in ^0[). One can also derive a one- dimensional integral representation (see 
eq. (26) of ||3^) which is valid for arbitrary e = (4 — n)/2 and therefore for arbitrary n, 



J(l,l,l) 



i vr^ ^ ?7 



1 f da a-'iiyay - (x/a)-') 



-piy+^ e 



[ya"^ -|- (1 — a; — y)a + x 



l-e 1 



(B.7) 



with X and y defined by eq. ( |B.6| ). 

Another way to get the result for arbitrary n is to use the connection between 
massless triangle integrals and the two-loop massive vacuum integrals, /(z/i, z/2, z/3). In 
particular, according to eq. (40) of [^, for arbitrary n the integral J(l, 1, 1) can be 
related (up to trivial leftovers) to the integral /(1, 1, 1), in such a way that one of the 
integrals is taken in 4 — 2e dimensions, and the other in 4 + 2e dimensions. Using this 



^^The main formula to be used is eq. (3.4) of ||2S 
are also presented in |2a|. 



Some explicit results for the integrals with Vi — 2 



29 



connection, and also known results for J(l, 1, 1) in arbitrary dimension |^2|, 0], one 
can reproduce the result for J(l, 1, 1). It can be written in terms of functions (see 
e.g. eq. (4.12) of 0). 

The result for J(l, 1, 1) in four dimensions is well-known |T3[ (see also in 
and can be presented in terms of dilogarithms, 



^{Pl,P2^P3) 

where 



71 



n=4 



^ I 2 [Lis (-px) +Li2 i-py)] + In ^ In + ln(pa;) \n{py) + y 

(B.8) 



1 — a; — ?/ + A 

or, in terms of the Clausen function, CI2 (see e.g. eq. (19) of similar representations 



are also given in |]42| , ^). In particular, in the symmetric case (see section 3.4) 



J(l,l,l) 



9 9 9 9 

Pi=P2=Pi=P 
n=4 



ivr (p{p ,p ,p ) 



i7r2 4 /TT 



n=4 



(B.IO) 



producing the same constant as the one (denoted as /) used in |Tl|]. The transcendental 
number CI2 (^^) = 1.0149417... corresponds to the maximum of Clausen's integral and 
appears frequently in two-loop calculations with masses (see e.g. [0, ^ ), the connection 
with the massless triangle diagrams being clear from ||41|| . 

If one is interested in expanding in e and calculating the integral J(l, 1, 1) up to order 
e, one can use eq. (29) of or eqs. (16) and (20) of |4l[| . 

In three dimensions, the result for the integral J(l, 1, 1) is very simple [^] and pro- 
-1/2 gg|^ ^Y^Q result around two dimensions, eq. (43) of can 



portional to {pfpip^ 

be used0. The corresponding two-dimensional integral has infrared singularities which 
can be regularized by the same e. However, the only non-trivial function is the same as 
in the four-dimensional case. The same equation (43) of ^l[] can also be used to get the 
results for higher values of n. 

When some of the external momenta squared vanish, the integral J(l, 1, 1) (considered 
in four dimensions) develops infrared (on-shell) singularities. The corresponding limits 
are considered in section 4. 



Appendix C: Two-point functions 

To check whether our results are consistent with the Ward-Slavnov-Taylor identity for 
the three-gluon vertex ( |2.12| ), we need expressions for the functions contributing to the 
gluon polarization operator and the ghost self energy. 

The corresponding scalar functions, J(p^) and G(p^), are defined in eqs. ( |2.7| ) and 
(|2.8|), respectively. The lowest-order results are J^^^ = G^^^ = 1. At the one-loop order, 
the results can be found e.g. in [Q. We present them here for completeness, and also to 

^^This is a special case of a more general result presented in . 
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show the proper normahzation of the functions. One-loop contributions to these functions 
are presented in Fig. 2. 

The gluon and ghost loop contributions to the J and G function in an arbitrary 
covariant gauge yield 



^2 



(47r)"/2 8(n -T) {^^^"^ ~ 2^ + ^^^"^ ' ^^^^'^ " " ^'^"^ - - 4)} 

(C.l) 



(C.2) 



while the quark loop contribution to the J function is 



^ ' 2NfTR 



{An 



|n/2 



n 



(C.3) 



(there is no quark contribution to the G function). The coefficients t], Ca and Tr are 
defined by eqs. ( |2.17D , ( |0| ) and ( |3.16| ), respectively. The function k(p^) is given by 
eq. i^). 

The ultraviolet divergences of eqs. ( |U.1|) -( |C75D are given by 



ja,uv) 



(An] 



2-e 



5 

3 ' 2'- 



(4vr; 



2-e 



Ca U2H)^ 



iuv) 



(C.4) 

where k^'^^-* = 1/e + . . . is the divergent part of k(p^), eq. ( |2.15| ). The corresponding 
counterterms are 



7t2 



Ca (I + |C 



q{1,CT) 



(4vr) 

(47r)2-^ 
t C 



1 



+ R)+0{e) 



(47r)2 4 



C, 



(4vr; 



2-e 



where g^ = g'^e~'^^{AnY , and i? is the renormalization-scheme constant chosen in such a 
way that i? = in the MS scheme (see also section 3.5). 

Note that in the Fried- Yennie gauge [00, ^ = —2, the ghost self-energy is finite 
as n — * 4. Moreover, if one chooses the ra-dimensional generalization of this gauge as 
^ = — 2/(n — 3) |5^, then the one-loop correction to the ghost self-energy vanishes. This 
is connected with the transversality of the gluon propagator (2^) in the coordinate space 
(at this value of ^). 

The gluon polarization operator is finite when 



71=4 



Ca 



(C.6) 



For Nf = 6, = i and Ca = 3, this value is = — ^ 
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This gauge was also used in refs. 
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Appendix D: Results for the ghost-gluon vertex 



There are two one-loop contributions to the ghost-gluon vertex which are shown in Fig. 3. 
Here, we present the most general results for the scalar functions contributing to the 
ghost-gluon vertex at the one-loop level. The definition of the ghost-gluon vertex and the 
decomposition in terms of scalar functions is given in eqs. ( p.9|) and (|2.10|) , respectively. 
The lowest-order expression is given by eq. (|2.11|) . In one-loop expressions, we use the 
notation for y), and rj which can be found in eqs. ( p.l4| ), ( p.l5| ) and ( p.l7| ), respec- 
tively. /C and Q denote the symmetric scalar combinations constructed from the external 
momenta, eqs. ( ^.2[ ) and (|3.1|). 

At the one-loop level (in the Feynman gauge) we get 



r(i?(pi,P2;P3) 



X 



+PVPl^3 [pI^PiP^)^ + {PiP2)i^I+pIi^2 + {P2P3)f^3 
-2P2^P2;.3 [pI{P2P3)^ + phi + {PiP2)k2 + (PlP3)/«3 
+Pl^P2^3 [- {pi{P2P3) + Plipm)) ^ + {PlP3)l^l + {P2P3)i^2+pIi^3 

+P2f,Pl^,, [{pliPm) + '2{P1P2){P2P3)) V + 3(PiP2)ki + 3pIk2 + 3{p2P3)l^3] } • (D-1) 

For the scalar functions (|2.1(]| ), eq. (p.l| ) yields 

a^^'°\p3,P2,Pi) = - ,^ 2/2 \ {(2p?-P2)v^ + /«i -2K2-K3} 



6(^'°)(p3,P2,Pl) = -^^^ ^ {p\{P2P3)^^ ^ P\l^\ + {P\P2)I^2 + {P\P3)l^3l , 



(D.2) 



(D.3) 



:(^'°)(p3,P2,Pi) 



r] 



r 1- 

(47r)"/2 ^ 4/C 



{{Pi - pI){PiP2)^ - {pI - {PiP2))ki + {pI - {piP2))i^2 + {pI - pI)i^3] , (D-4) 



d^'^'\p3,P2,Pl) = jl^,CA^ 



{(^2(^1^3) + 2(piP2)(P2P3)) + 3(piP2)/tl + 3pll^2 + KP2P3)i^3} , (D-5) 



X 

^2 



e^^'°\p3,P2,Pi) = - ^l^Y/2 ^ {^3(^1^2)^ + {PiP3)i^i + {P2P3)i^2+pIi^3} ■ (D-6) 

In the limit n A, the expressions for all scalar functions ( p.2| )-( ]IX^ ) have been 
compared with the corresponding results presented in [|l^ , taking into account the erratum 
(which affects the results for the c,d and e functions). The comparison was successful, 
with the exception of two minor things. One of them is related to the definition of 
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the renormahzation-scheme constant and was aheady mentioned before, see footnote in 
section 3.1. The second one is that in the erratum |13[ the sign of the term b{Pi, P2, P3) 
in the expression for the d function is changed from minus to plus (see p. 2554 of ||13|| ). 
However, the comparison is successful if we keep the original sign, which is minus^^. 

Furthermore, the results which follow from eq. ( |D.1| ) (contracted with pi) for two 
infrared-divergent cases, (i) Pi = pi = ^ and (ii) pi = P3 = 0, have been compared with 



those presented in table B.II of ref. |]T5|. The latter were obtained in the Feynman gauge 
and expanded in the limit n — > 4, keeping the finite (in e = (4 — n)/2) terms. To consider 
the limit of our expressions, we used formulae presented in section 4.3. We found that 
our results coincide in this limit with those from [jl5| . 



Now, we present the results for the scalar functions ( 2.10 ) in arbitrary covariant gauge: 



a(^'^)(p3,P2,Pi) 



(4vr)n/2 



16/C 



{ [2/C mn - 5) + 2) + e ((n - 3)^ + 2) pl{p,ps) - - 2)p^(piP2)] [(^2^3)^ + f^i 

+ [(e - 2) (4/c - ^pIpI) m + ^)pliP2P3)] [{pm)v + ^^2] 

+ [(^ - 2) ((^ + 2)/C - ^pI{p2P3)) m + 2) (P2P3)'] [{pm)^ + ^3] 

H{P2P3)^ [(e-2) {^ + 2pl{pm)) + ((n-3)e + 2) (/C-2(piP2)(P2P3))]} , 



(D.7) 



^ 'Hp3,P2,Pi) = -^^C^^^ 

((n - 4)e + 4) (n - 3)(pip2)(p2P3)/C 
+ {{n - 3)e + 2)p\pI (2(e - 2)/C - ^{n - l){piP2){p2P3)) 

Hii - 2)p2 ((n - 1)pIpI{p,P2) + 1C{pI - 2{n - 3)(piP2)))] [{P2P3)v + '^i] 
+ [-^ ((n-4)e + 4) (n - 3)(piP2)(P2P3)/C + e(^ - 2)pI {{n-l)plplpl + {n-2)m) 

+ ((n-3)e+2)p2 (]C{ipl+{?,^-A){pip2))-S,{n-l){p,p2f{p2P3))] [{pm)^ + ^^2] 
+ [-e ((n - 4)e + 4) (n - 4)(pip3)(p2P3)/C 

+ ((n - 3)e + 2)pI (e(P2P3)(/C - (n - l)(piP2)(piP3)) + 2{i - 2)K{pm)) 

- 2)pI {{n - l){p2P3){PiP3f - ^{pI + {P1P3)))] [{pm)'^ + 1^3] 
+^ [2i {{n - A)i + 4) (n - ^){piP2){piP3){p2P3)K: 

+ ((n-3)e + 2)p^ (2^(n-l)(pip2)'(piP3)(p2P3)-2/C(piP3)(2(^-2)(pip2)+eP2)-^^') 
H{i - 2)p2 (2(n - 2)lC{prP2){p2P3) - 2{n - \)p\pl{pW2){P2P3) + ^')] } > (D-8) 



(47r)"/2 16/C2p2 

x{ [((n-3)e+2) {i{n-l){p,p2f{pmf " /C (4(p?)2 + (piP3)((e + 2)p? + 3e(piP2)))) 



^^These misprints were confirmed by the authors of 
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+ 



+ 



- 2)(piP3) {{n - l){piP2){pmf + JC{2pl + 3(piP2)))] [{P2P3)<P + «i] 
((n - 3)e + 2) {apm){P2P3){{n - l){piP2? + 2)C) + 2X:p?(p2 - (p^p^))) 



+e(e-2)(piP3)(pf/C-(n-3)/CQ-( 
((n - 3)e + 2) - l)(piP2)(piP3)'(P2P3) 
-/C ((e - 2)p?(p? - p^) + e(p2P3)(2(piP3) - pD)) 

+^(^ - 2)(piP3)(p2P3) ((n - 3)/C - (n - I)p2p2 



l)PlP2P3) [(PiPa)^' + «2] 



[(PlP2)<^ + «3] 



+293 



((n - 3)e + 2) (P1P2) (2pf(pf - pi)IC + e(n - l){piP2){piP3r{p2P3) 
+C(P1P3)(P2P3) ((e-2)(n-l)(piP2)(piP3)' + e(^-2)/C(p2-p2)^]|^ 



(D.9) 



x{ ((^ - 4)C + 4) {piP2){piP3) {{n - 3){piP3){P2P3) + pI{PiP2)) 

- {{n - 3)C + 2)p3 {C{P2P3){pIpI - (n - 2)(pip2)(piP3)) + 6(pip2)/c) 

- 2)p3(p2P3)(Pl(P2P3) - (n - 2)(piP2)^)] [(P2P3)¥' + «l] 

+ - 4)^ + 4)(n - ?>){piP2){piP3){P2P3? 

+ ((n - 3)e + 2)p2 (e(piP2)(P2P3)((n - 3)(p2P3) - 2p^) - 4p2/c) 

-(e-2)p2 (e(p2P3)(p?P2 + («-l)P2(PlP2) + {Pmf) + (e-2)p^/C)] [(piP3)<^ + t^2] 

+ [e ((^ - 4)e + 4) (n - 4)(pip3)(p2P3)/C + e(C - 2)p^(piP2) ((n - l)p^(p2P3) - K) 
- ((n - 3)e + 2) p2(p2p3) (^(^ _ I)p2(p^p3) + + 6)/C)] [{pm)v + «3] 



+2(^ 



- ((n - 4)e + 4) (piP2)(P2P3)/C - 3)(piP3) + - 2)p^) 

+ ((n - 3)e + 2)p^ (e(piP2)(p2P3)(2/C + np^(pip3)) - '^{plf {pmf) 

+ (e - 2)p^(piP2)(piP3)(p2P3) {an - l)(piP2) + 2(e - l)pVj\ } , (D.IO) 



x{ [((n - 3)e + 2) (-e(n - l)p^(piP2)'(piP3) + /Cp?(2(pip3) - ^Pi 
- 2)(n - l)pl{piP2){piP3f + e{n - 2)JCpl (2pl + 3{pm 
((n - 3)e + 2) (P2P3) {-an - 1)pIpIpI + m(n - S)pl + 2pl)) 
+m - 2)pI {{n - 1)p\pIpI + n{n - 3)Q - pi) 
\{n - 3)C + 2) (-i{n - l)plplpl + /C((e + 2)p? _ - 3)^2) 



+ 



[(P2P3)¥' + 1^1] 



[{PIP3)V + «2] 



+ 



- 2)(P2P3) ((n - 1){P1P3? + 2/C 

+2(P2P3) [{{n - 3)e + 2) ipm) {^n - l)pl{PiP2f - 21Cpl) 



pI [iPlP2)^ + H3] 
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- 2)(n - l)plipiP2)iPiP3)' - - 2)}Cplipl - pj)] } .(D.ll 



From the expressions ( |D.7] )-( ]D.11| ), one can see that again, as in the three-gluon 
functions ( p.l0|) -( p.l5|) , the values ^ = and ^ = — 4/(ri — 4) are distinguished. Putting 
^ = 0, we get rid of the momenta squared in the denominators (only /C may remain, cf. 
eqs. ( p.2| )-( p^ ), while for ^ = —4/ (n — 4) one of the momenta squared, pi, still survives 
in the denominators of the c function, eq. ( p.9|) , and the e function, eq. ( p.ll|) . This does 
not matter, since (i) the c and e functions do not contribute to the Ward-Slavnov-Taylor 
identity for the three-gluon vertex, eq. ( p. 12 ), and (ii) for the proper ghost-gluon vertex, 
eq. ( p.9| ), we get an extra pf in the numerator as a result of contracting with pi^ (since the 
tensor structures corresponding to the c and e functions contain pi^, see eq. ( p.lO|) ). If 
one wants to put some of the momenta squared equal to zero in other gauges, one should 
carefully consider the appropriate limit and expand the functions in the numerator, in 
exactly the same way as it was described in Section 4 for the case of the three-gluon 
vertex. 

In the limit n — > 4, the only divergent function in the one-loop ghost-gluon vertex is 
the a function ( p.7| ), also since this is the only function which is present at the zero-loop 
level, eq. (|2.11| ). The ultraviolet-divergent part of the a function (|D.71 ) is 



a 



{1,UV) 



(UV) 



(D.12) 



(47r)2-= 

where, as in eq. (|3.37] ), k^'^'^^ = 1/6 + .. . corresponds to the divergent part of the function 
K, eq. (|2.15| ). Therefore, by analogy with (|3.38|) , the counterterm contribution is given 

byQ 



a 



(i,cr) 



Ca 

(47r)2 2 ' 



0[- + R 



{l-0[- + R]+O{e), (D.13) 



(47r)2-= 

where 'g^ = g'^ e"'^'^ {4:T^Y , and R is the renormalization-scheme constant (i? = in the MS 
scheme, which is achieved by a suitable extraction of the over-all factor in eq. ( p.l3|) ). In 
particular, there is no singularity in the one-loop ghost-gluon vertex in the Landau gauge 
(^ = 1), see e.g. ref. ||. 



Appendix E: Results for ^3 = in arbitrary gauge 

For arbitrary ^, the results for the gluon and ghost loop contributions to the scalar func- 
tions (^ areQ 

a^'^^\pIpI-.^) - " ^ 



(47r)"/2 64(n- l)(n-4) 

x|(n-l)(n-4)e(4+(n-4)0 

+ (n-4) [l6(3n-2) + 4^(n-l)(10n-33) - 5^^(n-l)(n-4)J [ki + K2] 

-lQ{n - 1) (n + (n - 3)0 (5i2)"' [ki - ^2]} , (E.l) 



pI 

Pi 



Pi 

P2 



^°See also the discussion of hdr, 9^, etc. in section 3.5. 

■^^We also consider non-equivalent permutations of the arguments. 
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(47r)'*/2 64(n - 1) {n - 4) (n - 6) 



Pi 



x ^ - l)(n - 4)(n - 6) e (4 + (n - 4)0 
\Pi 

+Ki{n - 4)(n - 6) [l6(3n - 2) + 4^(n - l)(10n - 33) - bi^{n - l)(n - 4)" 
+2/^2 (n - l)(n - 4)(n - 6) ^ (4(n - 3) - (n - 4)0 

-8(n - l)(n - 4) /«2 [((^i2)"'(n - 6) (2 + (n - 3)0 + ((^is)"' ^ + - 3)0 
-4(n-l) [«i-K2] [(5i2)"' (n-6)(4(n+l) + e(5n-18)-^'(n-3)) 

+ (5i2)-' (4 + (n - 4)0 ((n - 3)^ - n + 6) - 2 (^la)-'^ (n + (n - 3)0] } , (E.2) 

b^'^^\pIpIq)^-^^Ca ^ 



pI pI 

— o/^l o'^S 

Pi P2 



(47r)"/2 64(n- l)(n-4)(n-6) 
x|(n- l)(n-4)(n-6) C(4+ (n-4)0 

-8 (n - l)(n - 4) [^i + k^] {8^2)-^ ^{n + {n- 3)0 
+ [«i-«2] [l6(n-6)(4n2-21n+14) + 4^(n-l)(n-6)(10n^-77n+152) 
f (n - l)(n - 6)(5n' - 48n + 104) + 16 (^la)"' (n - 1)^ (n + (n - 3)0] } , (E.3) 



B''''\^.pIpI)--j^,Ca 



pI , o^'? 

-^Ki + 2^K2 

Pi P2 



(47r)"/2 " 64(n- l)(n-4)(n-6) 

x|-(n- l)(n-4)(n-6) C(4 + (n-4)0 

-8(n-l)(n-4)«2 [(n-6) (2 +(n-3)0 (l + ((^12)"') + ((^12)"' ^ (n+(n-3)0' 
-(n - 4)(n - 6) /ti [l6(4n - 3) + 20^^ - l)(2n - 7) - ^^'^{n - l)(n - 4)" 

+4(n - 1) [/ti - K2] [{5i2)-^ (n - 6) (4 + 3e(n - 2) + ^'(n - 3)) 
-((^12)-' (4 + (n - 4)0 ((n - 3)^ - n + 6) + 2(5i2)-' i{n+{n- 3)^)] } , (E.4) 



X ■ 



(47r)"/2 ^ 32(n - 1) (n - 4) (pf - pi) 
n-l)(n-4)e (4 + (n-4)0 
+ [ki-K2] [l6n(n-4)+40n-l)(6n^-41n+72)-3C^(n-l)(n-4)2]}, (E.5) 



pI pI 

-^K'l ^1^2 

Pi P2 



c^'^'\o,pi;pl) = 
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(47r)"/2 ^ 32(n - l)(n - 4)(n - 6) pf (pf - pi) 
X I 8(n - l)p? [ki - K2] [-{n - 6) (6 + (n - 5)0 



36 



+in - 6){5u)-' (6 + an - 7) - - 3)) - 2(5i2)-' Un + in- 3)0 
+8(n - l)(n - 4)(5i2)-' (p? + pI) K,^{n+{n- 3)0 



+ (n-4)(n-6)fi;i 



2\2 



(pi) 



Pi 



- 1)^ (4 + (n - 4)0 + 16 (np^ - (4n - 



+2e(n - 1) ((12n - 38)^^ - 4(5n - 17)p^) + f (n - l)(n - 4)(-3p2 + 4]?^) ] | , (E.6) 
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X • 



(47r)"/2 ^ 16(n - 1) (n - 4) {n - 6) {pi - pi) 

{n -l){n- 4)e (4 + (n - 4)0 (4^(71 - 3) - 5(n - 6)) 
2 



1 1 
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Pi P2 



2 2 
Pi -P2 



(n - 4) [ki + K2] - 2(^12) ^ [ki - K2] 



X [8(n - 6)(4n - 7) + 4^(n - l)(n - 6)(5n - 11) 

- l)(5n^ - 60n + 108) - 2^^n(n - l){n - 3)] } , (E.7) 



F^'^^\Q,plpl)--y^.CA 
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(47r)"/2 32(n - 1) (n - 4) (n - 6) p? (p? - pi) 

n-l)(n-4)e(4 + (n-4)0 



[n - 6)^/^1 + 4(n - 3)e^«:2 
Pi P2 



+{n - 4)ki [(32(n - 6)(3n + 1) + 4^(n - l)(n - 6)(20n - 63) 

-3C^(n - l)(n - 4)(n + 2) - 4^=^(n - l)(n - 3)(n - 4)) 
-{Su)~^ (32(n - 6)(n + 2) + A^n - l){7n + 30) 

- l)(3n2 + 50n - 24) - A^^{n - l){n - 3){n - 16)) 

-24(5i2)"^ (8(n - 6) + 8^(n - 1) + ^^(n - l)(3n - 2) + 2^3(n - l)(n - 3) 



+16p? 5 — 5 ^ [(^(^ - 6)(n2 - 12n + 17) - 2^(n -l){n- 6)(3n - 10) 



Pi - pi 



- l)(n - 5)(n2 - 6n + 12) + 2^^{n - l)(n - 3)(n - 5) 
+ {612)'^ (16(n - 4)(n - 6) + 4^(n - l)(4n - 15) 

+^'^{n - l)(n - 4)(5n - 3) + 3^^(n - l)(n - 3){n - 4)) 

+3(5i2)"' (8(n - 6) + 8e(n - 1) + f - l)(3n - 2) + 2C%n - l){n - 3))] } , (E.8) 



h^'''KpIpIo) 



5(2 T] 



(47r)"/2 



X' 





pi 


p? 


1- 


^Kl 


2 '^2 




Pi 


P2 



16(n — 1) (n — 4) (n — 6) (pf — pi) 

n - l)(n - 4)^ (4 + (n-4)0 (2^(^-3) - 3(n-6)) 

+2((5i2)"^ \{n - 4) + K2] - 2(5i2)"^ [/«i - K2] 
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X 



24(n - 6) + 2C{n -l){n + 12) + ^^(n - l)(lln - 12) + 6^^{n - l)(n - 3) 
+ [ki - K2] [lQ{n - 6)(n^-2n-2) + A^{n - l){n - Q){Qn^ -37n + m) 
-^\n - l){7n^ - 88n^ + 376n - 552) - 8^\n - l)(n - 3)(n - 5)] } , (E.9) 

where = n{pj), see eq. ( p.l5| ), while the coefficients r] and Ca are defined by eqs. ( p.l7| ) 
and ( |3.3|) , respectively. 



Appendix F: Results for pf = pl = in arbitrary gauge 

The scalar functions Ui corresponding to the decomposition of the three-gluon vertex in 
this limit are defined by eq. (|4.58| ) . Using the expressions for the one-loop contributions to 
the A, B, C, F and H functions, eqs. ( [4.43[ )-( P757[ ), and the representations ( [4.59| )-( ^^ ), 
we get the following results for the one-loop contributions to the Ui functions: 



(47r)"/2 ^ 8(n- l)(n-4) 
X {4{2n^ - 15n + 19) + 2^{n -l){n- 3)(4n - 17) - ^^{n -l){n- 4)^} , (F.l) 



1 



(47r)"/2 ^4(n-l)(n-4) 
X {2(3n2 - 17n + 8) + ^{n - l)(n - 4)(8n - 27) - ^^(n - l)(n - 4)^} , (F.2) 



K 



(47r)n/2 ^ 16(n- l)(n-4) 
X {4(n - 4)(n - 5) + 2^(n - l){2n'^ - 19n + 36) - ^^(n - l){n^ - 8n + 20)} , (F.3) 



(47r)"/2 4(n- l)(n-4)(n-6)p2 
X {6(n -2){n- 6)(n - 7) + ^{n - l){n - Q)(2n^ - 23n + 40) 

- l)(n - 2)(n2 - 15n + 57) - i^{n - l)(n - 3)(n - 8)} , (F.4) 

^ {4(n-4)2 + ^(n-l)(n-6)}, (F.5) 



(47r)"/2 ^ 4(r2- l)(n-4)p^ 



1 



(47r)"/2 16(n- l)(n-4)p2 
X {l6n(n - 4) + 4^(n - l)(6n^ - 41r2 + 72) - 3^^(n - l)(n - 4)^} , (F.6) 



^2 



1 



X 



(47r)"/2 ^ 16(r2- l)(n-4)p2 
{8(5n2-25n + 2) + 4^(n-l)(8n2-59n+112) -^2(n-l)(3n2-24n+40)}; (F.7) 
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g^V ^.rrp n{n - 3) 



(47r)"/2 



{n-l){n-2) 
n-2 



(4vr)"/2 



ANfTn 



n — 1 
1 



(4vr)"/2 
(4vr)n/2 / 



(n - l)(n - 2) 
n + 2 



(n -l){n- 2)p2 
n - 4 



5(2 



^2 



^ 4A^fTp 

(47r)"/2 ^ ^ (n - l)p 

71,2 — 4n + 6 



n-2 



2 ^' 



^2 r/ 



(47r)"/2 f ^ {n-l){n-2)f 



(F.8) 
(F.9) 
(F.IO) 
(F.ll) 
(F.12) 
(F.13) 
(F.14) 



where, as usual, k, = Comparison with the definition of the functions Fi given in 

eq. (29) of [jl4| shows that the functions f/i, f/2, t^s, ^4, f^s, Uq and f/7 are proportional to 
F2, F3, Fi, Fq, F5, F4 and F7, respectively. 
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Figure 1. One-loop three-gluon vertex diagrams. 




Figure 3. One-loop ghost-gluon vertex diagrams. 



